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ABSTRACT 

We present the first direct measurement of the angular-diameter distance with H2O mega- 
maser technique to a galaxy beyond 100 Mpc in a single step without any local calibration. We 
image the sub-parsec scale H2O maser disk at the center of NGC 6264 with four tracks of Very 
■ Long Baseline Interferometry (VLBI) observations. We presented an image from two of these 

tracks in MCP Paper III, and we add two new tracks here. With the 2.3 years of monthly Green 
Bank Telescope observations of the H2O emissions in NGC 6264, we measure the centripetal ac- 
celerations of the maser lines from the accretion disk by tracking the line-of-sight velocity of each 
maser component as they change with time. The measured accelerations suggest that the maser 
disk appears to have more substructure than the masers in UGC 3789 and NGC 4258. In order 
to model the more complex maser disk properly and determine a precise distance and Hubble 
constant, we adopt a Bayesian modeling program that fits the maser disk in three dimensions. 
By assuming the masers follow circular orbits around the central black hole, the best fit of the 
data gives an angular-diameter distance of 137±19 Mpc, a Hubble constant of Hq = 74 ±10 
km Mpc~^, and a mass of the central black hole of (2.89±0.39) x 10^ Mq. The result shown 
in this paper demonstrates the feasibility of the megamaser technique for measuring distances to 
galaxies deep in the Hubble flow. 

Subject headings: accretion, accretion disks - galaxies: nuclei - galaxies: masers - galaxies: active 
- galaxies: ISM - galaxies: Seyfert 



INTRODUCTION 



Cosmology research has entered a new era since the discovery of the acceleration of the Universe (Perl- 
mutter et al. 1999; Riess et al. 1998). "Dark Energy", which has negative pressure and accounts for 73% of 
the total energy density of the Universe, is currently the best candidate to explain the cosmic acceleration, 
and understanding its nature has become one of the most important problems in modern astronomy and 
astrophysics. 

There have been several very promising methods proposed to explore dark energy and its equation- 
of-state parameter w with high accuracy (Frieman, Turner, Huterer 2008). As pointed out by Hu (2005), 
among all observables for probing DE in light of the CMB, w is most sensitive to variations in _ffo- Hu (2005) 
concluded that the single most important complement to the CMB for measuring the DE cquation-of-state 



- 2 - 



parameter w at is a determination of the Hubble constant to better than a few percent. This insight 

forms the fundamental motivation for the Megamaser Cosmology project (MCP; Reid et al. 2009, Braatz et 
al. 2010), which aims to determine Hq to 3% accuracy. 

The key to a precise and direct determination of the Hubble constant is to measure accurate distances to 
galaxies well into the Hubble flow (i.e. > 50 Mpc). The galaxies must be distant to reduce the contribution 
of the uncertainty coming from peculiar velocities. Among all the approaches to measure precise distances 
directly, the megamaser method pioneered by the study of NGC 4258 (Herrnstein et al. 1999) has proven to 
be the most effective to make precise and direct distance measurements to galaxies beyond our Local Group. 
In the MCP, we bypass the extragalactic distance ladder and apply the megamaser technique to galaxies in 
the Hubble flow, obtaining direct angular-diameter distances in a single step without any local calibration. 
Thus, this approach enables us to obtain an independent measurement of the Hubble constant without 
appealing to the Extragalactic Distance Ladder. This was not possible in the past because of insufficient 
sensitivity to detect H2O megamasers at sufficiently large distances. In the past decade, the advent of the 
100-m Green Bank Telescope has made this project possible. 

The H2O megamaser method involves sub-milliarcsecond resolution imaging with Very Long Baseline 
Interferometry (VLBI) and single-dish monitoring of H2O maser emission from sub-parsec circumnuclear 
disks at the center of active galaxies. For a "simple" H2O maser disk such as the one at the center of NGC 
4258 (Herrnstein et al. 1999, Argon et al. 2007), which exhibits Keplerian rotation with the systemic maser^ 
mainly lying in a single narrow ring, one can in principle determine the geometric distance to the galaxy by 
measuring four orbital paramters of the H2O maser disk: 

D = -szn z , (1) 

ttsys 

where Yorb is the observed orbital velocity of the systemic masers, A9 and Qsys are the angular orbital radius 
and the observed acceleration of the systemic masers, respectively. In principle, Vo, and A9 can be inferred 
from the rotation curve of the accretion disk obtained from the VLBI imaging, i can be inferred with AO 
and the maser distribution in the sky (see Kuo et al. 2011), and agys can be measured with H2O maser 
spectra from multi-epoch monitoring with a single-dish telescope (see Braatz et al. 2010 for details). With 
the megamaser technique, Herrnstein et al. (1999) has measured a distance of 7.2 Mpc to NGC 4258 with 
7% uncertainty (including the systematic errors), and it is expected that, with the new VLBI data presented 
in Argon et al. (2007) and Humphreys et al. (2008), the uncertainty of the distance can be reduced to '^3% 
(GreenhiU 2009). 

In the Megamaser Cosmology Project, we aim to measure maser distances with ~10% or better accuracy 
to each of about 10 galaxies, thereby determining Hq to ~3% after averaging the results. The first galaxy 
measured by the MCP using the megamaser technique was UGC 3789, which has a distance of 49.9±7.0 Mpc 
(Braatz et al. 2010; Paper II hereafter). We expect to reduce the overall uncertainty to ~8% by including 
new data and better disk modeling (Reid et al. 2012). In this paper, we present the distance to NGC 6264, 
the first galaxy beyond 100 Mpc measured with the megamaser method. In section 2, we present our VLBI 
and single-dish observations. In section 3, we show the analysis of the centripetal accelerations of the masers 
in NGC 6264. The analyses of the Hubble constant and distance determination are presented in section 
4. In section 5, we discuss the chanllcnges of applying the maser technique to distant galaxies. Finally, we 
summarize the results in section 6. 



^the maser components having velocities close to the systemic velocity of the galaxy. 
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2. Observations and Data Reduction 
2.1. VLBI Data 

We observed NGC 6264 with four tracks of VLBI observations made between 2009 and 2010 using 
the Very Long Baseline Array (vlbaH augmented by the 100-m Green Bank Telescope (GBT) and the 
Effelsberg 100-m telescopeH. Two of the four tracks (BB261F and BB261H) were presented in Kuo et al. 
(2011), and in this paper we present the observations BB261K and BB261Q. Our final VLBI image includes 
the data from all four tracks. The observing parameters for the last two VLBI tracks are shown in Table 
1. 

The observation and data reduction procedures for BB261K and BB261Q are essentially the same as 
described in section 2.2 in Kuo et al. (2011). Here, we only discuss the slightly different approach to 
self-calibration used for the two observations presented in this paper. 

Typically, in order to achieve sufficient sensitivity in the self-calibration process, we require a single 
maser line (linewidth~2 km s~^) with flux density above ^80 mJy, or a sufficient number of maser lines in 
a single intermediate frequency (IF) band with flux densities above ~30 mJy. In the latter case, we achieve 
high enough signal-to-noise ratio (SNR) by averaging multiple maser lines in narrow ranges of both velocity 
and space to perform phase calibration. However, in BB261K and BB261Q, the maser lines were generally 
weaker than those in the previous observations on NGC 6264, and even averaging multiple maser lines in a 
single IF band was not sufficient to achieve sufficient SNR for self-calibration. 

One way to overcome this difficulty is to utilize all maser lines with flux densities above ^20 mJy in all 
IF bands for self-calibration. To do this, we flrst remove the phase shifts due to different maser locations 
by dividing a model of maser positions into the data. The model division effectively moves all maser spots 
in consideration to the phase center of the observation, and since these masers now have the same phase, 
averaging the maser lines in the calibration process will increase the SNR directly. 

To choose the model, the best approach would be to use as the initial guesses the maser positions 
obtained from phase-referencing observations of the same source, if available, or from observations at certain 



^The VLBA is a facility of the National Radio Astronomy Observatory, which is operated by the Associated Universities, 
Inc. under a cooperative agreement with the National Science Foundation (NSF). 

^The Effelsberg 100-m telescope is a facility of the Max-Planck-Institut fiir Radioastronomie 



Table 1. Observing Parameters 



Experiment 
Code 


Date 


Galaxy 


Antennas'^ 


Synthesized Beam 
(mas X mas,deg)^ 


Sensitivity 
(mJy)"^ 


Observing 
Mode" 


BB261K 
BB261Q 


2009 Nov 25 

2010 Jan 15 


NGC 6264 
NGC 6264 


VLBA, GB, EB, VLA 
VLBA, GB, EB 


1.10x0. 30, -9. 7 
0.80x0.24,-11.8 


0.5 
0.5 


Self-cal. 
Self-cal. 



""VLBA: Very Long Baseline Array; GB: The Green Bank Teleseope of NRAO; EB: Max-Planck-Institut fiir Radioas- 
tronomie 100 m antenna in Effelsberg, Germany; VLA: The Very Large Array of NRAO. 

*^The average FWHM beam size and position angle (PA; measured east of north) at the frequency of systemic masers. 

"^"Self-cal." means that the observation was conducted in the "self-calibration" mode of observation. 

'^Tlie sensitivities were measured with a channel width of 3.6 km s~^. 
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cpochs when the self-cahbration can be done easily because of stronger niascr fluxes. Following this guideline, 
for BB261K and BB261Q, we used the averaged maser positions from BB261F and BB261H, and divided the 
model into the data in order to shift each maser spot to the phase center before performing self-calibration 
observation. Once the model subtraction and line-averaging are done, we treat the masers as a single maser 
spot, and follow the same calibration procedure for solving the atmospheric phase, solution editing, and 
imaging as decribed in Kuo et al. 2011. Finally, we combine the maser map obtained from BB261K and 
BB261Q with the one from BB261F and BB261H on the image plane by aligning the maser disks using the 
average position of the systemic maser features from BB261K and BB261Q with velocities between 10191.1 
km s~^ and 10198.3 km s~^ as the reference. We present the combined data in Table 2. 



2.2. GBT monitoring 

We took 20 H2O maser spectra with the GBT for NGC 6264 between 2008 November 21 and 2011 
March 2. Except during the summer months when the humidity makes observations at 22 GHz inefBcient, 

we took a spectrum on a monthly timescale. For these observations, we follow the same observing settings 
and data reduction procedures as in Braatz et al. (2010). Table 1 shows the observing date and sensitivity 
for each observation. A representative H2O maser spectrum for NGC 6264 can be found in Figure 1 in Kuo 
et al. 2011. 
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Table 2. Sample data for NGC 6264 



V * 

» op 




(56,^ 


B,b 








(km s-i) 


(mas) 


(mas) 


(mas) 


(mas) 


(km yi:~^) 


(km s""'^ yr""'^) 


1 r»m o o o 
iUyiO. 00 


A OAT 

U.o97 


A AA C£; 

U.0056 


A AO O 

-O.Ooo 


A A 1 /I 1 

0.0141 


A AA 
0.00 


1 A A 
i.OO 


iUyi4. / 1 


U.4Uc5 


A AA^C 


A A1 1 
U.Ull 


A A1 C /I 

U.Ulo4 


A AA 
U.UU 


1 AA 
l.UU 


1 noi 1 no 
iUyii.Uy 


r\ A AT 
U.4U ( 


A AA/I Q 

U.UU40 


A AAC 

-U.UUo 


A A1 1 A 


A 1 1 
U.ll 


1 AA 
i.UU 


1 AAA7 /I T 

iUyU7.47 


A O AA 

U.399 


A AA/?A 

U.0069 


A AOf? 
0.U2D 


A A 1 Tl 
O.Ui ( i 


A AA 
0.00 


1 A A 
i.OO 


1 AO C O 1 C 

iUooo.io 


A A OT 


A AAOT 

U.00o7 


A A 1 1 

-O.Uil 


A A A A 1 
0.0091 


A AA 
0.00 


1 A A 
i.OO 


1 AO /I n CO 

11)849.00 


A Ad A 

U.494 


A AAOi; 

0.0036 


A AAO 

O.UUo 


A AAOA 
0.0089 


A AA 
0.00 


1 AA 
1.00 


1 AC/1 p: m 


A r^Ai 
U.oUi 


A AA01 


A A07 


A AACC 


A 1 
U.iO 


1 AA 
i.UU 


1 AO /I o on 


A ^A/1 

U.0U4 


A AA1 


A AOQ 


A A A /in 
U.UU4y 


A A /I 

-U.U4 


1 AA 
i.UU 


1 A O O O £^ ^7 


A C AO 


A AA /I O 
U.0042 


A AO A 

-0.U24 


A Alio 

O.Ollo 


^1 O /I 

-0.24 


1 A A 
i.UU 


1 Aoo /I 1 n 


A COT 


A AATI 
U.007i 


A AOA 

-o.uoy 


A A1 CA 

0.0159 


A 

0.27 


1 AA 
l.UU 


10809.70 


0.550 


0.0037 


-0.028 


0.0092 


0.00 


1.00 


10806.08 


0.561 


0.0034 


-0.053 


0.0090 


0.00 


1.00 


10802.46 


0.572 


0.0017 


-0.031 


0.0044 


-0.09 


1.00 


10798.84 


0.577 


0.0013 


-0.040 


0.0034 


-0.27 


1.00 


10793.40 


0.582 


0.0014 


-0.028 


0.0039 


-0.05 


1.00 


10787.97 


0.596 


0.0046 


-0.032 


0.0125 


0.00 


1.00 


10780.73 


0.614 


0.0031 


-0.045 


0.0075 


-0.13 


1.00 


10773.49 


0.632 


0.0058 


-0.025 


0.0149 


0.00 


1.00 


10766.25 


0.654 


0.0051 


-0.046 


0.0136 


-0.26 


1.00 


10762.63 


0.673 


0.0055 


-0.088 


0.0142 


0.00 


1.00 


10755.38 


0.664 


0.0044 


-0.071 


0.0117 


-0.06 


1.00 


10751.76 


0.672 


0.0020 


-0.052 


0.0052 


-0.07 


1.00 


10748.14 


0.677 


0.0035 


-0.084 


0.0099 


0.02 


1.00 



Note. — Sample of data for NGC 6264. The entirety of data for all galaxies is 
available in the electronic version. 

'^Velocity referenced to the LSR and using the optical definition (no relativistic 

corrections). 

''East-west and north-south position offsets and uncertainties measured relative to 
the average position of the systemic masers from the BB261K and BB261Q data. 
Position uncertainties refiect fitted random errors only. 

'^Measured or estimated acceleration and its uncertainty for each maser component. 
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Table 3. GBT Observing dates and sensitivities for NGC 6264 



Ep.uli 




Date 


Dav AiiiiilKT 


'l'..,„. (iv) 


riiis 


PcTicKl 





2008 


November 21 





30.9 


1.7 


A 


1 


2009 


January 16 


56 


28.2 


1.3 


A 


2 


2009 


February 3 


74 


35.9 


1.3 


A 


3 


2009 


March 4 


103 


31.9 


1.2 


A 


4 


2009 


March 31 


130 


37.8 


1.4 


A 


5 


2009 


May 13 


172 


41.5 


1.7 


A 


6 


2009 


November 7 


351 


44.3 


1.7 


B 


7 


2009 


December 12 


386 


27.4 


1.5 


B 


8 


2010 


January 11 


416 


26.3 


1.1 


B 


9 


2010 


February 9 


445 


39.0 


1.5 


B 


10 


2010 


March 7 


471 


28.3 


1.5 


B 


11 


2010 


April 13 


507 


36.8 


2.2 


B 


12 


2010 


May 10 


534 


28.2 


1.2 


B 


13 


2010 


July 2 


588 


40.3 


2.2 


B 


14 


2010 


October 30 


708 


40.2 


1.7 


C 


15 


2010 


November 26 


735 


38.4 


1.5 


C 


16 


2010 


December 24 


763 


30.3 


1.5 


C 


17 


2011 


January 23 


793 


31.4 


1.2 


C 


18 


2011 


February 7 


808 


42.7 


1.4 


C 


19 


2011 


March 2 


831 


39.0 


2.3 


C 



Note. — The rms noise values arc calculated without performing Manning smoothing to the 
spectra and are based on 0.33 km s~^ channels. We label Period A, B, and C to those times 
when we have continuous observations on a monthly timescale. These periods are separated 
by summer months during which the humidity makes observations inefficient. 
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3. Acceleration Analysis 
3.1. Methodology 

The accelerations of the H2O maser lines measured in the megamaser technique are centripetal accelera- 
tions from circular orbits in the accretion disks. We measure these accelerations by tracking the line-of-sight 
velocity of each maser line as they change with time. We adopted two different methods to determine the 
accelerations of maser spots in this paper. The first method is what we call the Eye-tracking method. In 
this method, we first plot the velocities of maser lines as a function of time (e.g. the left panel in Fig. 1), 
and identify a group of data points within a particular time segment by eye that we believe to be from a 
particular maser line that persists over time. We then fit the trend of data points in the plot with straight 
lines. The slopes of the fitted lines directly give the accelerations. 

The second method is the GLOFIT method. Here, we use a global least-square fitting program (Braatz 
et al. 2010; Humphreys et al. 2008) to fit multiple maser lines from multiple epochs simultaneously, and 
determine the accelerations of the fitted maser lines in a single step. This is the same method we used for 
measuring the maser accelerations in UGC 3789 (Braatz et al. 2010; Reid et al. 2012) 

For NGC 6264, we adopt the GLOFIT method to measure accelerations of the systemic masers and the 

Eye-tracking method for the high velocity mascrs. Wc choose the GLOFIT technique for the systemic masers 
because accurate acceleration measurements of the systemic masers is critical for a precise Hubble constant 
determination with the megamaser technique, and it is the GLOFIT method that can better handle the line- 
blending effect of the systemic masers and determine their accelerations more accurately. On the other hand, 
since many high velocity maser lines arc relatively isolated, and very precise acceleration measurements of 
these masers are not critical for accurate distance determination, to avoid the time-consuming and laborious 
fitting process for all high velocity maser lines, we apply the Eye-tracking method to provide good estimates 
of accelerations for the high velocity masers. 

While it is more appropriate to adopt the GLOFIT program to measure the accelerations of the systemic 
masers, it may be necessary to use a different strategy to apply the method in NGC 6264 than in UGC 3789 
because as a result of the joint effect of more significant blending and SNR limitations, finding stable and 
converging solutions is more difficult. One way to improve the situation is to fit all the maser lines in the 
velocity range of interest to drift with a single acceleration in the global fitting. In principle, this approach 
should help the program to reach convergence more efficiently and stably because the SNR of the fitting 
increases when we enforce a common acceleration. We have found that this is indeed the case after testing 
this approach with synthetic spectra generated from simulations. 

Finally, the applicability of the GLOFIT program in the way described above depends critically on 

the assumption that masers within the selected range of spectra drift at the same rate. Therefore, to 
avoid systematic errors it is important to examine whether the chosen sections of spectra really satisfy this 
assumption. A correct verification of the constant acceleration assumption would depend on many factors 
including whether the masers show the same linear trend in the position-velocity (P-V) diagram (indicating 
they come from the same radius), the persistency of the maser spectral pattern with time, the reduced 
of the fit, and the degree of line- blending. For the interested readers, we refer to Kuo 2012 (Ph.D. thesis) 
for details. 
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3.2. Acceleration Fitting for High Velocity Meisers 

In Figure 1, we plot the radial velocities of NGC 6264 maser peaks measured by eye as a function of 
time. For the blueshifted and redshifted masers, we first identify lines that are persistent in time and then fit 
a straight line to the data to measure the acceleration directly with the eye-tracking method. We estimate 
the uncertainty of the measurements by scaling the fitting error by the square root of reduced x^- In Table 
4, we show the measured accelerations and uncertainties. 
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^*t^ — ■ 





+ ++ + 
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+ +++■ 



+ + 



1 1 



4*J 



400 
Doy 



Fig. 1. — In this figure, we plot the radial velocities of NGC 6264 maser peaks as a function of time (the 
crosses) along with the fitting results from the eye-balling method (for the high velocity masers) or from the 
modified GLOFIT program (for the systemic masers). The data between Day and 200 come from spectra 
taken in Period A; the data between Day 300 to 600 from spectra in Period B; and the data between Day 
700 to 900 are from spectra in Period C. 



The weighted average accelerations of the redshifted and blueshifted masers are -0.06 and 0.01 km s~^ yr~^ 
respectively. The rms scatter of acceleration of the redshifted and blueshifted masers are both 0.11 km s~^ yr~^ 
The small acceleration and rms scatter indicate that the high velocity masers are close to the mid-line of the 
accretion disk as expected. 

When assigning the measured accelerations to the corresponding channels in the VLBI datascts. we pay 
attention to whether the periods over which the accelerations are measured actually bracket the epochs of 
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Table 4. Acceleration Measurements for the High Velocity Masers in NGC 6264 



Velocity 
(km s~^) 


Acceleration 
(km yr~^) 


(km yr~^) 




n no 
— u.uy 


u.u / 


O/ion AO 




u.oo 




nit; 

U. AO 


n /ion 


yoo'i.yo 


Q 


n 1 7 
u. 1 / 


yoou.oy 


n 1 Q 
u. ly 


12 


yoyo. 10 


n 1 Q 
— u. lo 


U.UfS 


you 1.04 


n Oft 


U.UfS 


Qfli 1 '7r\ 
yoii. ( u 


n 1 1 
u. 1 1 


n 1 n 

U. lU 


Qft1 A '7C\ 


n 1 
— u. 1^ 


U.UD 


yoou.yy 




U.Uo 


yooo.^o 


— U.U4 


u.u / 


yoo ( .00 


n no 
u.uy 


n 1 n 

U. lU 


yo / y.^y 


u.^y 


U.Uo 


9692.04 


0.17 


0.10 


9787.36 


0.02 


0.03 


9811.77 


0.05 


0.06 


10749.05 


0.02 


0.14 


10751.24 


-0.07 


0.04 


10753.28 


-0.03 


0.08 


10755.20 


-0.10 


0.10 


10766.50 


-0.26 


0.03 


10776.50 


-0.30 


0.42 


10781.77 


-0.13 


0.02 


10793.18 


-0.05 


0.10 


10798.96 


-0.27 


0.11 


10802.19 


-0.09 


0.08 


10809.58 


-0.17 


0.06 


10824.78 


0.27 


0.34 


10826.46 


-0.75 


0.73 


10837.53 


-0.16 


0.08 


10840.21 


-0.27 


0.06 


10841.58 


-0.04 


0.05 


10845.09 


0.16 


0.06 


10848.68 


0.00 


0.11 


10873.30 


0.00 


0.16 


10912.41 


0.11 


0.60 



Note. — The acceleration measurements for the 
high velocity mascrs. The components having ve- 
locity higher than 10745 km s "'" are redshifted 
masers whereas those having velocity lower than 
9815 km s~^ are blueshifted masers 
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our VLBI observations. We have two tracks (BB261F & BB261H) of VLBI observations in Period A (2008 
November 21 - 2009 May 13; see Table 2), two tracks (BB261K & BB261Q) in Period B (2009 November 
7 - 2010 July 2) and no data in Period C (2010 October 30 - 2011 May 2). For those masers persisting 
over Periods A through C, we assign the measured accelerations to the corresponding VLBI channels of the 
data from averaging all four tracks. For those masers that persist only during Period A or B, we assign 
the accelerations only to the VLBI channels from observations in the respective period. For high velocity 
maser spots without acceleration measurements, we assume their accelerations are zero and assign them to 
the data from averaging all four tracks. 

3.3. Acceleration Fitting for Systemic Masers 

Figure 2 shows the maser spectra from Epochs through 5 (Period A) and from Epochs 6 through 11 
(Period B). We do not show the spectra from Period C because there is evidence that the systemic masers 
in this period have accelerations different from masers in Period A and B. Since we do not have VLBI data 
for Period C, the acceleration measurements for the systemic masers in this period arc not useful for the 
distance determination, and thus we do not show the analysis here. Since both Periods A and B bracket two 
tracks of VLBI observations and the accelerations of systemic masers could change with time, we measure 
acceleration in Periods A and B separately. In cases where the accelerations measured in both periods within 
a certain velocity range are the same and the corresponding VLBI positions are consistent, we average the 
accelerations and VLBI positions. For those maser lines that only persist during one period, we only use the 
VLBI positions at this particular period for these masers in the distance and Hubble constant measurement. 

To measure the accelerations, we break the spectra into distinct velocity ranges, as indicated in Figure 
2, and analyze each range individually. 

3.3.1. Clump 1 

The results from the eye-tracking fitting suggest that the maser lines at V ~ 10203 km s^^ have 
significantly smaller acceleration (a ^ 0.6 km s~^ yr^^) than the rest of the masers (a ~ 1.2 km s~^ yr^^). 
Therefore, to fit a constant acceleration for each velocity range, we divide Clump 1 into two velocity sections: 
10183.0 to 10200.2 km s"! and 10200.2 to 10203.5 km s'S where V = 10200.2 km s"^ is at the local minimum 

that divides the two sub-clumps. 

For the maser lines of Clump 1 within 10183.0 to 10200.2 km s^^, we could not find a model to make a 
stable fit unless we narrow the velocity range to 10183.0 to 10197.7 km s^^. This suggests that maser lines 
between 10197.7 and 10200.2 km s~^ have different accelerations than the other lines in the fitted velocity 
range. Since we cannot give a good constraint on the accelerations of these lines, we do not include them in 
the distance determination. For lines between 10183.0 and 10197.7 km s~^, we fit 8 Gaussian components to 
the spectra and obtain an acceleration of 1.04± 0.14 km s~^ yr~^. The reduced for this fit is 1.144 (153 the 
degrees of freedom). For the maser lines between 10200.2 to 10203.5 km s~^, we fit two Gaussian components 
and obtain an acceleration of 0.75±0.10 km s^^ yr^^. We overplot the acceleration measurements in this and 
the following subsections on top of the time- velocity plot for the systemic masers (the middle plot of Figure 
1) and list the results in Table 4. We also show an example of the result of the Gaussian decomposition in 
the acceleration fitting with the masers in Clump 1 in Figure 3. 




Fig. 2. — The upper panel shows the spectra from epoch through 5 (purple, blue, green, yellow, and 
orange), and the bottom panel shows the spectra from epoch 6 through 11 (purple, blue, green, yellow, and 
orange). The velocity ranges of the systemic masers in Period A and B are divided into 7 clumps for the 
convenience of acceleration measurement. 



3.3.2. Clump 2 

The maser lines in Clump 2 are apparently drifting and the general pattern of the clump seems to 
persist with time. However, we couldn't find a good model to fit the data. The best reduced we can 
achieve is 1.44, and the fitted acceleration is 1.2±0.3 km yr~^- The situation does not improve when 
we only fit subsets of this clump. The large reduced suggests that our constant acceleration assumption 
may not hold. Therefore, for masers in this clump, we only use the VLBI channel at V=10218.11 km s^^ 
in the distance determination because it is the only line for which we have an acceleration measurement 
(a=1.55±0.46 km s~^ yr~^) from the eye-tracking fitting for this clump. 



3.3.3. Clump 3 

Maser lines in Clump 3 have the largest acceleration among the systemic masers in NGC 6264. The 
acceleration is large enough that the drifting of the whole clump can be seen clearly in Figure 2. We fit ten 
Gaussian components with an average linewidth of 2.4 km s~^ to the masers between 10226.0 to 10252.3 
km s~^. The measured acceleration is 4.43±0.36 km s^^ yr~^ with a reduced of 1.006 (314 degrees of 
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freedom). 

3.34. Clump 4 

The maser lines in Clump 4 reside within the same velocity range as Clump 1, but cover time period 
B, and the situation is similar. So we perform the fitting in two velocity sections: 10184 to 10201.6 km 
and 10201.6 to 10205.5 km s~^, where V = 10201.6 km is at the local minimum that divides the two 
sub-clumps. For mascr lines between 10184 to 10201.6 km s~^, wc fit 9 Gaussian components to the data 
and obtain an acceleration of 1.12±0.18 km s^^ yr^^ with reduced of 1-009. 

For the masers between 10201.6 to 10205.5 km s~^, we could not find a reliable fit because of severe 
blending. By comparing the maser spectra from Period A and Period B, we argue that since the line structure 
remains similar over Periods A and B at this velocity range, wc arc seeing the same maser lines that appear 
at slightly different velocities because of the acceleration. This idea is supported by the fact that the VLBI 
position of the masers in this velocity range from Periods A and B are very well consistent with each other. 
For these reasons, we assume these masers have the same acceleration (a = 0.75±0.1 km s~^ yr"^) as the 
masers in Clump 1 at the similar velocity range. 

3.3.5. Clump 5 

Clump 5 covers a velocity range between 10207 to 10212 km s~^ and the lines within this velocity range 
do not show a clear sign of drifting. In addition, the line structure changes substantially over the course of 
time and we found it difficult to measure a reliable acceleration for this clump. We can measure significantly 
different accelerations, ranging from 0.4 to 4.2 km s~^ yr^^ by using different subsets of the spectra in the 
fit. This implies that we may see different lines at different epochs. Along with the fact that the VLBI 
positions for masers in this clump have relatively large error bars, we decide not to include these masers in 
the Hubble constant determination. 



3.3.6. Clump 6 

The acceleration measurement for masers in clump 6 is also not trivial, like Clump 5. We couldn't 
achieve a converging fit. The line structure changes substantially at epoch 10. By inspecting the spectra in 
Period B along with those in Period C, wc discovered that the line at V ~ 10224.1 km s"'^ starts to drift 
at a much higher rate (a > 7 km s~^ yr^^) after epoch 9. This may be caused by newly arising lines from 
smaller radii of the disk. Because of this complexity, which could introduce significant systematic error in 
our fitting caused by maser lines with different accelerations, we do not include this group of masers in the 
Hubble constant determination. 
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Model and Data 




101S4 101S6 10IS8 10190 10192 10194 10196 
Velocity (km/s) 



Fig. 3. — An example of the Gaussian decomposition for the acceleration measurement. In this example, we 
fit the masers between 10183.0 and 10197.7 km in the spectra. The panels from top to bottom show the 
spectra (black lines) from epoch through 5. Each of the eight Gaussian components fitted to the data are 
represented by different colors. The purple curves at the bottom of each panel are the residuals from the fit. 



3.3.7. Clump 7 

The masers in Clump 7 cover nearly the same velocity range as Clump 3, but the drifting of the whole 
maser pattern becomes much less clear because of both severe blending and reduced SNR. The blending 
causes the line structure to become very smooth and there are only two lines that can be seen to drift clearly 
in the time-velocity plot (Figure 2). For this clump, we fit nine Gaussian components with an average 
linewidth of 2.3 km s""'^ to the masers between 10230.5 to 10250.5 km s~^. The measured acceleration is 
3.96±0.59 km s^^ yr~^ with a reduced of 1.075 (180 degrees of freedom). This is consistent with the 
value for clump 3, and implies that masers in clumps 3 and 7 are basically the same maser emission but 
appear at different times. 
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4. Distance Determination for NGC 6264 and the Hubble Constant 

4.1. The 3-Dimensional Bayesian Approach to Hubble Constant Determination 

To determine the Hubble constant with the NGC 6264 data, we adopt a Bayesian approach to model 
the maser disk in three dimensions (Reid et al. 2012). Unlike the "ring method" used in previous result for 
UGC 3789 (Braatz et al. 2010), the Bayesian technique can fit a maser disk with maser clouds distributed 
arbitrarily across the thin disk. 

The 3-dimcnsional nature of the modeling allows one to model the warping structure and orbital incli- 
nations precisely if there is a sufficient number of both high velocity and systemic maser spots. Note that 
accurate modeling of the disk structure not only helps to measure the impact parameters (i.e. the projected 
separation between a maser spot and the BH) of maser spots more precisely, but also is crucial to place tight 
constraints on the position of the black hole and the recession velocity of the galaxy, both of which are very 
helpful for constraining the maser distance and Hubble constant. 

The Bayesian nature of the analysis allows one to use all prior information for the model parameters 
to constrain the fitting more easily. In principle, the use of prior information can help guide the fitting 
to find the most likely solutions efficiently and reliably in the large space which could have multiple 
local minimums. Furthermore, rather than seeking a single (best) solution and formal uncertainties, the 
Bayesian approach can directly probe the posteriori probability distribution of parameters without the 
assumption that the uncertainties follow a Gaussian distribution. Finally, one can model the eccentricity of 
the maser orbits in three dimensions, and directly measure its impact in the distance determination with 
the H2O megamaser method. However, since our data is currently too noisy to allow a convergent fit for 
the eccentricity parameters, in this paper we will only model the maser disk with circular orbits. We will 
explore the effects of eccentricity in a future paper after data quality improves by including new data in the 
analysis. 

The fitting program uses the same code for fitting the maser disk in NGC 4258 (Humphreys et al. in 
prep) and UGC 3789 (Reid et al. 2012) In this program, we fit the maser data including the East & North 
offset (x & y) of the maser positions on the sky, and the observed radial velocities (v) and accelerations (a) 
of the maser spots with a three dimensional model of the maser disk. We construct the model (Figure 4) by 
placing each maser spot at radius r and disk azimuth angle (f) (with cj) for systemic masers and (j) ~ 90° 



Tabk^ 5. Accc4c^ration ]\I(^asur(;ments for thc^ Systc^mic Masers in NGC ()2G4 



Clump 


Velocity Range 


Epochs 


Num. of 


Lincwidth 


Acceleration 




xl 


d.o.f. 




(km s"^) 




Components 


(km s^^) 


(km s~^ yf"^) 


(km s~^ yr"^) 






1 


10183.0 - 10197.7 


0-5 


8 


1.9 


1.04 


0.40 


1.144 


153 


1 


10200.2 - 10203.5 


0-5 


2 


1.9 


0.75 


0.14 


0.966 


40 


2 


10209.1 - 10228.0 


1-5 


11 


1.7 


1.20 


0.30 


1.440 


144 


3 


10226.0 - 10252.3 


0-5 


10 


2.4 


4.43 


1.13 


1.006 


314 


4 


10184.0 - 10201.6 


6-11 


9 


2.0 


1.12 


0.54 


1.009 


183 


7 


10230.5 - 10250.5 


7-11 


9 


2.3 


3.96 


1.77 


1.075 


180 



Note. — Col(l): The clump number; Col(2) The velocity range for acceleration measurement; Col(3) The epochs of the spectra 
used for fitting; Col{4) The number of Gaussian components that fit the data; Col(5) The average linewdith of fitted lines; Col(6) 
The best fit acceleration; Col(7) The uncertainty of the acceleration; Col(8) The reduced of the fit; Col(9) The number of degrees 
of freedom of the fit. 
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for redshifted masers) on an elliptical orbit about a point mass (presumably a black hole) Mbh, with the 
whole system at a distance D. The center of the maser disk on the sky is {xq, yo), its recession velocity is 
Vgys, and the velocity correction to the pure Hubble flow is Ucor- 

The disk can be warped in two dimensions: the inclination warp i(r) and the position-angle warp p{r), 
each of which is specified by 3 parameters: 

iir) = ia + kr + i2r^ , (2) 

P{r) ^ Po + Pir + P2r'^ , (3) 

where r is the radius of the orbit in milli-arcsec. One can choose how many of these warping parameters one 
should actually use in the fitting depending on the degree of warping. During the fitting process, the program 
adjusts r and along with the above-mentioned global parameters including Hq and Mbh simultaneously 
until the best fit is reached. Note that in choosing the set of parameters to model the maser disk, unlike 
what was done for modeling the maser disk in NGC 4258 in Herrnstein et al. (1999), we fit Hq instead of 
using D as one of the primary parameters. Details are given in the appendix. 

The program adopts a Markov chain Monte Carlo (McMC) approach (e.g. Geyer 1992; Gilks, Richardson 
& Spiegelhalter 1996) to obtain the posteriori probability distribution function of the model parameters. In 
this approach, one fits the data in a "brute force" manner by randomly trying a large number of model 
parameter values. Because of the "Markov chain" nature of the parameters, the z*'* -|- 1 trial parameter 
values are dependent on, and close to, the i*'' values. In particular, this program uses the Metropolis-Hastings 
algorithm (e.g. Metropolis et al. 1953; Hastings 1970) to choose the Markov-chain trial parameters, and 
a global step-size factor in the fitting is adjusted to scale parameter steps so that an optimal Metropolis- 
Hastings acceptance rate of w23% (Gelman, Gilks, & Roberts 1997) is achieved. This algorithm has the 
property that the distribution of the trial parameter values is equal to the desired posteriori probability 
distribution function of the parameters. In the actual fitting process, we first run 20 burn-in stages, each 
with 10^ trials, in order to arrive at near optimum parameter values. After discarding the burn-in stage 
trials, we evaluated 10^ McMC trials to obtain final posteriori probability distributions. Note that the trial 
number is choosen to be such a large number in order to ensure we explore the parameter space sufficiently. 



4.2. Hubble Constant Derived Prom the Bayesiann Fitting 

The top panel of Figure 5 shows that the degrees of both the inclination warp and the position-angle 
warp appear to be relatively small from the maser map. Therefore, we fit the warp only to the linear term in 
Equations (3) & (4). To account for systematic uncertainties in our measurements of the positions, velocities, 
and accelerations of the maser spots, we estimate more realistic "error floors" of 8 mas and 16 mas for x & y, 
respectively, 0.8 km s~^ and 0.4 km s~^ for velocities of systemic and high- velocity maser features, and 1.0 
km s~^ yr^^ for accelerations of high-velocity masers. For accelerations of systemic masers, the estimated 
error floors range from 0.3 to 0.7 km s~^ yr^^, depending on each individual maser clump. We add the error 
floors to these observational errors in quadrature before fitting the data. 

We adopt flat priors for all parameters except for iJp (72±18 km s^^ Mpc~^), Vgys (10196±28 km s""'^; 
Beers et al. 1995)0, and Vcor (— 19±300 km s ^). The first two priors were purposely set fairly "loosely" so 



*A11 velocities quote here are non-relativistic, optical definition, LSR velocities; these velocities were converted internally in 
the fitting program to rclativistically correct values. 
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as not to overly constrain and dominate the fitting. The prior for Vcor comes from the sum of two values: a 
correction of —19 km s~^ km s~^ to convert from the LSR to the cosmic microwave background frame, and 
an assumed value of the peculiar velocity of the galaxy expected from large scale gravitational perturbations. 
Since we do not have an estimate of the peculiar velocity from simulation of large scale cosmic flows, we 
assume the peculiar velocity to be km s~^ with 300 km s~^ uncertainty. 

Among the priors adopted for the Baycsian fitting, the Hq prior especially plays an important role to 
get our final result. By cxporing the fitting with a flat prior on Hq and weighting the data differently, we 
discovered that there are two primary local minima in the space for the Hq parameter: one at Hq ~74 
km s~^ Mpc~^ and the other at Hq ~90 km s~^ Mpc~^. Without prior information on Hq, it is difficult 
to tell from the data itself which local minimum has a higher probability, and the fitting result can depend 
sensitively on how we weight the data, even though the weights may only change slightly. Introducing a 
prior on Hq allows the fit to converge to the local minimum of Hq '~74 km s~^ Mpc~^ more stably with little 
infiuence by changes in the data weights. We expect that this situation will improve in the future when we 
incorporate new data for this galaxy. 

Figure 5 shows the qualities of the Baycsian fitting by comparing the best fit model with the observed 
maser distribution, position- velocity diagram, and acceleration measurements. In general, the model matches 
the observations well, and the differences between the data and model are consistent with real uncertainties. 
The total reduced of fit is 0.951. 

The left panel in Figure 6 shows the maser distribution from overhead. As expected, the systemic masers 
reside at different radii because of their divergent accelerations, and the high velocity masers arc located 
close to the mid-line of the maser disk because their accelerations are close to zero. The latter feature 
substantially supports the assumption made in Kuo et al. (2011) that the high velocity masers are close 
to the mid-line of the disk within ~13 degrees, and therefore the BH mass can be reliably and accurately 
constrained by fitting a Keplerian rotation curve to the high velocity masers. In the right panel of Figure 6, 
we show the warp structure with the maser spots plotted on top of the warping disk. Note that we increase 
the disk inclination by 5 degrees in order to show the warp structure more clearly. 

We show the posteriori probability distribution for Hq from the Bayesian fitting in the left panel of 
Figure 7. The posteriori probability distribution is slightly asymmetric about the peak of the distribution, 

which has a 68% confidence range of of 9 km s""'^ Mpc""'^ centered at Hq =74 km Mpc""'^. The mode 
and the median of the posteriori probability distribution are 73 km Mpc~^ and 74 km s~^ Mpc~^, 
respectively. 

We found that the best fit Hq has some sensitivity on the error fioors assigned for the measured velocities 
of the masers. Currently, we adopt an error ffoor of 0.4 km s~^ for the velocities of high- velocity masers 

and 0.8 km for those of systemic masers. These numbers are estimates based on our understanding and 
experience of the effects of astrophysical noise and degree of maser line blending, and could be uncertain for 
about a fewxO.l km s~^. By exploring the variations in the best-fit Hq with different error fioors for the 
maser velocities, the fitted Hq could change as much as 5 km s^^ Mpc^^. We consider this as systematic error. 
After including this value in the total uncertainty budget, oiu- data gives a Hubble constant of Hq =74±10 
km s~^ Mpc~^. We show the rest of the fitting results in Table 5. 
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4.3. The Angular-diameter Distance to NGC 6264 

In the modeUng of the maser disk, we use Hq instead of D as one of the fitting parameters, as explained 
in the appendix. In order to obtain the posteriori probability distribution for the maser distance; to NGC 
6264, we calculate D with the equation D = Vsys/Ho for each Markov-chain trial in the Bayesian modeling, 
and compute the posteriori probability distribution for D afterward. The right panel of Figure 7 shows 
the posteriori probability distribution for D. The distribution appears to be asymmetric with a long tail 
toward the large distance end. This tail is a natural result of Hq being in the denominator of D = Vsys/Ho, 
with Hq having a substantial and roughly equal 68% uncertainty on either side of the peak of the posteriori 
probability distribution. The posteriori probability distribution of D has a 68% confidence range of 16 Mpc 
centered at D =137 Mpc. Since Hq has a systematic error of ~5 km s~^ Mpc~^ caused by the error floor 
uncertainties for maser velocities, we add the corresponing systematic uncertainty in the distance (9 Mpc) 
in quadrature to the error derived from the posteriori probability distribution, and the total uncertainty 
becomes 19 Mpc. 

5. Discussion 

5.1. The First Geometric Distemce to A Galaxy Beyond 100 Mpc 

Traditionally, Hubble constant determinations relied on indirect distance measurements through the 
Extragalacitc Distance Ladder (EDL), which bases its calibration on the Period-Luminosity (PL) relation 
of Cepheid variables. In the past, the need to appeal to the EDL was inevitable because direct distance 
determinations were limited to stars in the Milky Way, the Large Magellanic Cloud, and water masers in 
NGC 4258 (Distance = 7.2±0.5 Mpc; Herrnstein et al. 1999). Direct distance measurements to galaxies in 
the Hubble flow (Distance > 50 Mpc) have never been achieved until recently (e.g. Braatz et al. 2010; Suyu 
et al. 2010; Reid et al. 2012). 

In this paper, we obtain the first direct measurement of angular- diameter distance to a galaxy beyond 
100 Mpc in a single step without any local calibration. NGC 6264 is more than 15 times farther than NGC 
4258, and the angular size of the maser disk is more than 10 times smaller. In addition, the maser flux 
densities are over 100 times fainter. The small size (^ 1.4 mas in diameter) and the low flux densities of 
these maser disks have not only been a big issue for efficient VLBI imaging and accurate astrometry, but also 
make reliable acceleration measurements difiicult because of low signal-to-noise ratio. Moreover, instead of 
residing in a single narrow ring as in NGC 4258, the systemic masers in NGC 6264 are located in multiple 
radii. This prevents using the simple approach adopted for measuring the distance to UGC 3789 (Braatz 
et al. 2010), and instead we use the more sophisticated 3-dimensional modeling for this megamaser disk. 
The work presented in the previous MCP papers (Braatz et al. 2010; Kuo et al. 2011), Reid et al. (2012), 
and this paper shows the effort to face these challenges by developing the 3-dimensional Bayesian modeling 
program and the new ways to perform VLBI self-calibration and global acceleration fitting. It is because 
of these efforts that make sensitive VLBI imaging and accurate distance/Hubble constant determination 
possible. 



6. Summary 

Our main conclusions are the following: 
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1. Wc modeled the H2O maser disk in NGC 6264 with Bayesian analysis, and obtain a Hubble of Hq = 
73.3 ±10.3 km s~^ Mpc~^ and an angular-diameter distance of 136.5±16.3 Mpc. The Hubble constant 
is consistent with the value obtained in Paper II {Hq = 69± 11) as well as the Cepheid-based Hq from 
Preedman & Madore (2010; Ho = 73± 4) and Riess et al. (2011; Ho = 74± 2). This is the first time 
the megamaser technique has been successfully applied to a galaxy beyond 100 Mpc. 

2. The application of the megamaser technique to galaxies deep in the Hubble flow is intrinsically more 
difficult because of significantly lower flux densities and smaller disk angular size, and these factors 
make precise astrometry and acceleration measurements more challenging. In order to image the maser 
disks in distant galaxies more efficiently, utilizing multiple maser lines at different locations of the maser 
disk with the self-calibration technique developed in this work to enhance the signal-to-noise ratio may 
often be necessary. 

3. The measured accelerations suggest that the maser disk in NGC 6264 appc^ars to have more substructure 
than the masers in UGC 3789 and NGC 4258. The systemic masers are located at least at five different 
radii from the BH. Full 3-Dimensional modeling is essential to measure the distance to NGC 6264 
because of the complexity of the maser disks. 

In the future, we expect that with the new observations included in the analysis, we would be able to 
improve the accuracy of the distance measurement to better than 10%, and explore the impact of orbital 
eccentricity on the distance measurement 

A. Correlation Between Distance or Hq and Black Hole Mass Estimates 

In our Bayesian modeling of maser disks in the MCP, we fit Hq directly, rather than fitting D as in 
Herrnstein et al. (1999). There are four reasons that make Hq preferable over D in the disk modeling: 

(1) since Hq is the primary parameter we intend to measure, it is more convenient to save an extra step of 

calculation and have the posteriori probability distribution of Hq directly output from the modeling program; 

(2) in this way the effect of uncertainty in peculiar velocity of a galaxy can be directly incorporated in the 
Hq calculation; (3) since D oc l/ttgys and Ugys usually has the largest fractional uncertainty among our 
observables, as a result, the posteriori probability distribution function of D would usually be asymmetric 
even when the data quality is high; finally, (4) if we use D as a fitting parameter, its posteriori probability 
distribution is usually highly correlated with Mbh (See Figure 8a), typically about 0.99 or even 1. 

Such a high correlation coefficient between D and Mbh implies that these parameters may be degenerate, 
and the fitting result in such a circumstance might not be reliable. This problem was not recognized in 
previous work on NGC 4258 by Herrnstein et al. (1999). We found that replacing D with Hq in Bayesian 
modeling resolves this issue; the absolute value of the correlation coeffecient between Hq and Mbh is 
significantly less than one. Nonetheless, the nearly perfect correlation between D and Mbh deserves some 
further investigation in order to understand its nature and whether the results from the work on NGC 4258 
and from our initial modelings that fit the maser disks with D and Mbh are reliable. In the following 
analysis, we use a "simple" maser disk with systemic masers lying in a single narrow ring to demonstrate 
that the nearly perfect correlation between D and Mbh is not a sign of degeneracy, and is just a result of 
substantial discrepancies between uncertainties of different observables. 

In the problem of modeling the maser disk and determining the angular-diameter distance, the rela- 
tionship between D and Mbh arises from two different measurements : the centripetal acceleration of the 
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systemic masers and the Keplerian rotation curve for the high velocity masers. In the former case, we can 
see the relation between D and Mbh in the equation for the centripetal acceleration Usys- 

_ GMbh _ GMbh f.^. 

"•^VS — ^2 ~ £)2^02 ' ^■^^) 

where G is the gravitational constant, r is the physical radius of the orbit of the systemic masers, and A9 is 
the apparent angular radius of the systemic masers. In principle, can be obtained from the intersection 
points between the straight line traced by systemic masers {V = fl 9) and the Keplerian rotation curve traced 
by the high velocity masers (V = k ^~^/^) in the position-velocity diagram of the maser disk (see Braatz et 
al. 2010). By solving for with the rotation curves for the systemic and high- velocity masers, we obtain 
A6 = [k/Q)'^/^. With this, we can re-write Equation (Al) and express Mbh as a quadratic function of D : 

M,„ = 2gi(i)"'D^ (A2) 



On the other hand, since the coefficient k of the Keplerian rotation curve is defined as k= GMbh /D, 
we can also find a linear relationship between D and Mbh as 

Mbh = D. (A3) 

In Figure 8b, we plot Equation (A2) with black quadratic curves and Equation (A3) with blue straight 
lines for a given set of measurements of agys, k, and fl including their uncertainties. Every intersection of 
a blue line and a black curve (red dots in Figure 8b) represents a solution for D and Mbh from fitting the 
maser disk, and all possible solutions (within one a uncertainty) only fall in a region on the D — Mbh diagram 
whose shape is defined by the relative sizes of the imccrtainties in asys, k, and fl. In current megamaser 
galaxies being observed for distance determination, the uncertainty in k is usually negligible (~1%) while the 
errors in agys and SI are more significant (e.g. >10%). As a result, the blue lines span a very narrow region 
on the D — Mbh diagram while the black curves cover a much broader region, and all the possible solutions 
for D and Mbh fall on a region whose shape resembles a line. This makes D and Mbh appear to be nearly 
100% correlated. Therefore, the apparent high correlation between D and Mbh from the modeling program 
is not a result of intrinsic degeneracy, and thus the fitting results should be reliable. 

In Figure 8c, we demonstrate how the high correlation between fitting parameters is removed by replacing 
D with Vsys/Ho in Equations (A2) & (A3). Now, one can see that the possible solutions no longer fall exactly 
on a straight line, and actually follow the Mbh oc 1/Hq relationship. Because of this, the nearly perfect 
correlation doesn't exist between Mbh and Hq. Nonetheless, one can imagine that if the uncertainties in 
Usys and O are also small, then the black curves may cover a much smaller region, and the possible solutions 
may only fall in a small segment of the 1/Hq curve. The distribution would thus resemble a straight line, 
and high parameter correlation appears again. This is indeed seen in our data, but because of the analysis 
shown here, high correlation should no longer be a serious concern. 
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3D Bayesian Fitting 
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Fig. 4. — The model and parameters for Bayesian fitting of maser disks in the MCP. The cartoon pfot on 
the left shows the model maser disk around a supermassive BH at the center. The model maser spots are 
allowed to be at any locations on the disk plane. The fitting program can fit elhptical orbits, but in this 
paper we model the masers in circular orbits. In the plot, r and denote the radius and disk azimuth 
angle of a maser spot with cj) oi 0° for systemic masers (green spots) and (p ~ -90° or 90° for rcdshifted (red 
spots) or blueshifted (blue spots) masers. The 15 parameters on the right are the parameters that one can 
use to fit a maser disk. Among these parameters, Hq-, Mbh^ Vsys, xq, yo, i, and p are the most important 
and necessary for Hubble constant determination. Whether other parmameters should be included in the 
modeling would depend on degree of disk warping, circularity of maser orbits, and whether one wants to 
include uncertainty of galaxy peculiar velocity in the calculation. 
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NGC 6264 H20 Masers 
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Fig. 5. — Data (colored dots) and best-fit model (lines and black dots). Top panel: Positions on the sky. 
Middle panel: LSR velocity versus position along the disk. Bottom panel: Accelerations versus impact 
parameter. 
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Model Maser Distribution The Warped Maser Disk 
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Fig. 6. — The left panel shows the model maser distribution in NGC 6264 from the overhead perspective. 
The right panel shows the best-fit warp from the observer's perspective with model maser spots plotted on 
top of it. For illustration, we decrease the disk inclination by ^5° to show the degree of disk warping more 
clearly. 




Hubble Constant (km/s/Mpc) Distance (Mpc) 

Fig. 7. — The left and right panel show the posteriori probability distribution of the Hubble constant Hq 
and the distance to NGC 6264, D, respectively. The posteriori probability distribution of Hq has a 68% 
confidence range that centers at Hq =74.3 km Mpc~^ with an uncertainty of 9.0 km s~^ Mpc^^ while 
that of D has a 68% confidence range that centers at D =137 Mpc with an uncertainty of 16 Mpc. 
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Table 6. NGC 6264 H2O Maser Model 



Parameter 


Priors 


Posterioris 


Units 




72±18 


74±10 


km Mpc~^ 


^sys 


10196±28 


10209±2 


km s~ 


VcOT 


-19±300 


-24+291 


km s~^ 


Mbh 


,'5.45 


2.89±0.39 


10'' M© 


xo 


0.006 


0.008±0.003 


mas 


yo 


0.006 


0.007±0.003 


mas 


io 


88.5 


88.8±2 


deg 


ii 


1.9 


1.0±2.7 


dcg mas""*" 


Pn 


8."?. .3 


84.0 -■- 1.5 


dog 


I'L 


IS. 8 


IT.S . 2.1 


d,.gma^ ' 



Note. — Parameters are as follows: Hubble constant {Hq 
)) Vlsr of central black hole (Vsj^s), peculiar velocity with 
respect to Hubble flow in cosmic microwave background frame 
(V'cor)) black hole mass (Mbh), eastward (xq) and northward 
(yo) position of black hole with respect to the average position 
of the systemic maser features from the BB261K and BB261Q 
data, disk inclination (io) and inclination warping (change of 
inclination with radius: ii = di/dr ), disk position angle (po) 
and position angle warping (change of position position angle 
with radius: pi = dp/dr ). Flat priors are listed without 
uncertainties; otherwise Gaussian sigmas are listed. 
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Fig. 8. — We illustrate the apparent correlation between D and Mbh and the possible cause. The left panel 
shows the fitting results for galaxy distance D and BH mass Mbh from our initial version of the Bayesian 
fitting program plotted on the D — Mbh diagram. The nearly straight-line appearance shows that D and 
Mbh are highly correlated, and the correlation coefficient is usually about 0.99 or even 1.0. The plot shown 
in the middle panel provides an explanation for this correlation. The black curves represent the quadratic 
relation between D and Mbh from Equation (A2) and the blue straight lines show the linear relation for D 
and Mbh from Equation (A3) for a given set of measurements of Osys, and fi. For a given distance, the 
black curves (from top to bottom) give the Mbh measured with agys, k, and Q, plus its 1 ct, cr and -1 a 
uncertainty. Likewise, the blue lines (from top to bottom) give the measured Mbh based on Equation (A3) 
plus its 1 tJ, (7, and -1 a uncertainty. Every intersection of a blue line and a black curve (the red spots) 
represents a solution for D and Mbh from fitting the maser disk. From the discussion in the appendix, 
we conclude that it is because the uncertainty in k is usually significantly smaller than Ugys and il that all 
solutions fall nearly on a straight line on the D — Mbh diagram. The presentation in the right panel is 
similar to that of the middle panel, with the only difference that D in Equation (A2) and (A3) is replaced 
by Vsys/Ho, and we plot the Hubble constant in the x-axis. It can be seen that the solutions no longer fall 
on a straight line, and Hq and Mbh arc not perfectly correlated. 



